Co-existing chiral and collinear phases in a distorted triangular antiferromagnet 
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The entire magnetic phase diagram of the quasi two dimensional (2D) magnet on a distorted 
triangular lattice KFe(Mo04)2 is outlined by means of magnetization, specific heat, and neutron 
diffraction measurements. It is found that the spin network breaks down into two almost independent 
magnetic subsystems. One subsystem is a collinear antiferromagnet that shows a simple spin-flop 
behavior in applied fields. The other is a helimagnet that instead goes through a series of exotic 
commensurate-incommensurate phase transformations. In the various phases one observes either 
true 3D order or quasi- 2D order. The experimental findings are compared to theoretical predictions 
found in literature. 



Helical magnetic structures in frustrated spin networks 
have been known for several decades, see, e.g, |l|, |2|, [si, l4j . 
The interest in this phenomenon was recently rekindled 
by a discovery of a new class of multiferroic compounds, 
where ferroelectricity is inherently connected to chiral 
incommensurate magnetic order (e.g. la,^). The latter 
enables a coupling between the ferroelectric order param- 
eter and a magnetic field, usually forbidden due to their 
respective symmetries. Perhaps the simplest, yet funda- 
mentally the most important model that can provide the 
necessary chiral state is the two-dimensional triangular- 
lattice antiferromagnet (TLAF) that has a well-known 
"120°" spin structure. 

A spectacular example of a multiferroic behavior due to 
the TLAF environment was recently found in the layered 
molybdenate RbFe(Mo04)2 [7|. In that material trian- 
gular planes of almost classical Fe^+ S = 5/2 spins are 
well separated by non-magnetic M0O4 layers. The Fe'^''" 
subsystem is an almost perfect realization of the TLAF 
model [8|, Isl]. In RbFe(Mo04)2 an external magnetic field 
disrupts the chiral magnetic structure and thereby re- 
moves the multiferroic effect. Even more interesting be- 
havior can be expected in materials with a slightly dis- 
torted TLAF spin network. The distortion will unbal- 
ance the 120° structure and produce an incommensurate 
planar spin spiral. In an external magnetic field, theory 
predicts a multitude of exotic incommensurate and com- 
mensurate, chiral and centrosymmetric phases |2|. The 
chiral states could, at least in principle, possess multifer- 
roic properties. The focus of this work is KFe(Mo04)2, 
a material structurally similar to RbFe(Mo04)2 but one 
that actually realizes the distorted TLAF model. Us- 
ing a combination of experimental techniques, we find 
a remarkably complex magnetic phase diagram with co- 



existing collinear and helical structures, ordered in either 
three or two dimensions. 

At high temperatures KFe(Mo04)2 has a perfect trian- 
gular spin network due to the crystal symmetry of D^^ 
group with a = 5.66 A and c = 7.12 A. The distor- 
tion occurs as a result of a crystallographic phase transi- 
tion at T = 311 K fid, llll|. The low-temperature phase 
is monoclinic, with a doubling of the period along the 
c axis. Each Fe-plane becomes a distorted triangular 
lattice with two unequal exchange constants Ji and J2 
(Fig. [1]). Moreover, the adjacent Fe-layers become crys- 
tallographically inequivalent, with exchange constants J( 
and Jj. The actual lattice distortion is too small to be 
detected with the resolution of our experiments, and we 
shall henceforth adopt a hexagonal lattice notation. A 
previous ESR study [I4I of KFe(Mo04)2 came to a seem- 
ingly paradoxical conclusion: at low temperatures, in 
zero field, a helical spin structure coexists with a collinear 
state. It was hypothesized that the two type of magnetic 
order reside almost independently in the two inequivalent 
types of distorted triangular spin lattices, referred to as 
"S-layers" and "C-layers" , respectively. 

To verify this bold assumption we performed magnetic 
neutron diffraction experiments using single-crystal sam- 
ples from the same batch [12|. The crystals are trans- 



parent thin plates, of a natural triangular shape, typi- 
cally 1-5 mg, with the planes of natural growth perpen- 
dicular to the 3-fold axis. The data were taken on the 
HB-1 and HB-IA 3-axis spectrometers at ORNL operat- 
ing in 2-axis mode, using a Pyrolitic Graphite PG(002) 
monochromator to select A — 2.46 A for HB-l and 
A — 2.37 A for HB-lA. At low temperatures, two sets of 
magnetic Bragg peaks emerge, with propagation vectors 
(1/3 - C,l/3 - CO), C = 0.038, and (1/2,0,0), respec- 
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FIG. 1: (a) Arrangement of magnetic ions in the crystal struc- 
ture of KFe(Mo04)2 showing two inequivalent Fe^^ planes, 
(b) Schematic representation of the zero-field magnetic struc- 
ture. 



tively. The observed ordering temperatures for the two 
sets of reflections are identical within experimental accu- 
racy: Tn = 2.4 K. An analysis of 19 inequivalent Bragg 
reflections in the {h,k,0) plane at T = 1.5 K revealed 
that the (1/3 — C,, 1/3 — C,,0) peaks can be entirely ac- 
counted for by a planar helimagnetic state, with spins 
rotating in the (a, b) plane. At the same time, the 14 in- 
equivalent sets of (l/2,0,0)-type Bragg intensities mea- 
sured in the {h, fc, 0) plane are consistent with a coUinear 
AF spin arrangement, with spins in the (a, b) plane and 
forming a small angle of 15° with the a axis. Thus the 
diffraction data confirm the original two-layer model de- 
picted in Fig. [T|d. As discussed in Ref. [I4], the drasti- 
cally different spin arrangement in C- and S-layers can 
be accounted by the difference in the corresponding ra- 
tios R — J1/J2 vs. R' — J'll J2- theory predicts a switch 
from a helimagnetic to a coUinear state at i? > 2 [21. 

The complexity of the H — T phase diagram of 
KFe(Mo04)2 was revealed in bulk magnetic and calori- 
metric measurements. Steady-state magnetization data 
were collected in fields up to 12 T using an Oxford In- 
struments vibrating sample magnetometer and specific 
heat data were collected on a Quantum Design PPMS at 
Warwick University. High-field magnetization data were 
taken in the fields up to 25 T using a pulsed magnet 
at the KYOKUGEN center. Typical experimental C{T) 
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FIG. 2: (Color online) (a) Field derivative of the magnetiza- 
tion measured at T = 1.6 K. Solid arrows are experimentally 
observed anomalies. Dashed arrows and spin diagrams are 
explained in the text, (b) Specific heat measured on a Img 
KFe(Mo04)2 single crystal sample in magnetic fields applied 
in the (a, b) plane. 



and dM/dH curves are shown in Fig. [2l In zero field we 
observe a sharp specific heat anomaly at Ti = 2.5 K. In 
a magnetic field applied in the (a, b) plane this anomaly 
shifts to lower temperatures and survives up to 7? = 5 T. 
For H > 2 T an additional peak is observed at T2 > Ti. 
Beyond H > 5 T the Ti anomaly is replaced by a new 
feature at T3. The latter also moves to lower temper- 
atures with increasing field. Yet another maximum in 
specific heat is observed at T4 < T2 in the high field 
regime. In the magnetization data, for a field applied 
along the c axis, the only observed feature is saturation 
£^t -ffsatll = 16.9 T. However, for a field in the (a, b) plane, 
the magnetization curves show four distinct anomalies. A 
jump of the field derivative of magnetization dM/dH at 
Hi ^ 1.2 T is followed by sharp maxima at H2 — 4.5 T 
and i/3 ~ 7.5 T. Finally, a saturation is reached at 
Hsa,t± = 14.5 T. For different orientations of the field 
within the (a, b) plane, and for samples with varying pop- 
ulations of the three crystallographic domain types, the 
measured magnetization curves are very similar. In par- 
ticular, the angular variation of H2 and H3 is 10% and 
7%, respectively. Some samples show a smeared peak in 
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FIG. 3: (Color online) Cumulative magnetic phase diagram of 
KFe(Mo04)2 for a magnetic field applied in the {a,b) plane. 
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FIG. 4: Main panels: field dependencies of magnetic Bragg 
intensities measured in KFe(Mo04)2 at T = 100 mK. Insets: 
typical /-scans measured across the corresponding reflections. 



dM/dH in the field range above 8 T in the temperature 
interval 2 K< T < 2.5 K. This anomaly, presented by the 
magnetic field H^ in Fig. [31 was not visible in the sample 
used for pulse measurements. The magnetic and thermo- 
dynamic anomalies described above, together with the 
M{T) data from Ref. 12], allow us to reconstruct the 
entire H — T phase diagram, as shown in Fig. [3l The 
low-field part of the phase boundary at T2 is uncertain 
because of the absence of this anomaly at H < 2 T. 

The propagation vectors in the various phases were 
determined in neutron experiments on the D23 lifting 
counter diffractometer at ILL using a PG(002) monochro- 
mator and A = 2.38 A neutrons. The field was ap- 
plied along the b axis (in hexagonal notations). Sample 
environment was a dilution refrigerator. Typical mea- 
sured field dependencies of Bragg intensities are shown 
in Fig. m At T = 100 mK the intensities of commen- 
surate (1/2,0, 0)-type peaks go through two consecutive 
jumps a,t Hi = 1.3 T and iJ{ = 2.1 T. These transi- 
tions seem to have no effect on the incommensurate re- 
fiections. At T = 1.5 K a single intensity .jump is de- 
tected at if 1 = 1 T. As previously discussed in [I4I, the 
transition is to be associate with a spin flop in the C- 
planes, all moments rotating to be perpendicular to the 



applied field. The additional transition seen at low tem- 
perature at H[ requires further investigation. At higher 
fields all the action occurs within the S-planes. The in- 
tensity of the incommensurate (1/3 — C, 1/3 — Ci 0)-type 
refiections decreases and vanishes beyond H2 ~ 4 T. 
Within experimental resolution the value of the magnetic 
propagation vector is field- independent. Beyond H2 the 
(1/3 — C, 1/3 — C; 0)-type peaks are replaced by commen- 
surate reflections of type (1/3, 1/3, 0). The latter first in- 
creases in intensity, peaks at around 6 T, and decreases 
at higher fields to vanish at H3 ^^ 8 — 9 T. At still higher 
fields, no magnetic reflections were found on either the 
{h, h, 0), (1/3 - C, 1/3 - C, 0: or (1/3, 1/3, reciprocal- 
space rods. At the temperature of neutron measurements 
(T — 100 mK and T = 1.5 K) we found no signature of 
a high-field transition that could be associated with the 
T4-anomaly described above. 

A remarkable feature of the neutron data is the dif- 
ferent dimensionality of magnetic ordering in the dif- 
ferent phases. Scans across the (1/2, 0,0)-type and 
(1/3 — ^7 1/3 "Cj 0)-type refiections are resolution-limited 
along the h-, k- and /- directions. In contrast, the 
(1/3, l/3,0)-type peaks in the regime H2 < H < H3 
are actually Bragg rods parallel to the c axis, stretching 
across much of the Brillouin zone in the /-directions. The 



corresponding c-axis correlation length is only 6 lattice 
units. The measured correlation length within the (a, b) 
plane is much larger, about 100 lattice units. The 2D 
character of the ordering persists in samples cooled in a 
6 T applied field. The transition to the short range or- 
dered state at T2 naturally demonstrates a smeared C{T) 
anomaly in contrast to a sharp peak at the transition to 
the 3D ordered phase at Ti. 

Guidance to understanding the complex phases real- 
ized in the 5'-_planes can be drawn from the theoretical 
work of Ref. [3|. At first, we can estimate the relevant 
exchange parameters from the saturation fields and sus- 
ceptibility. Because the neutron reflections observed in 
high fields correspond to C-planes, we assume the sat- 
uration is associated with this kind of planes, though 
S-planes, naturally, also give a contribution to the mag- 
netization. For the C-planes, the saturation fields are 
given by: 5/iB-ffsati = 8(J{ + J2)S and 5/iB-ffsat|| = 
8( J{ -|- J2)S + 2DS, respectively, where D is the single- 
ion easy-plane anisotropy defined as in [9|]. Using the 
observed values of Hsat± and i?sat|| we get J{ -I- J2 = 0.96 
K and D' = 0.32 K. Beyond the spin-flop at Hi, the 
C-planes contribution to magnetic susceptibility should 
be constant: xc = 3^Mb/[8(>^i + >^2)]- Subtracting this 
value from dM/dH data for Hd < H < Hc2 we ob- 
tain the susceptibility of the S-planes: xs — 0.12 /^b/T 
per Fe^+ ion. Using the latter value and the measured 
C = 0.038, by applying the equations in Ref. [3|, we get 
Ji — 0.37 K; J2 ~ 0.69 K. As a self-consistency check, 
the experimental ratio R — 0.53 warrants a helimagnetic 
ground state for the S-layers in zero field [4I . 

Now, the measured values C and Hsati. can be applied 
to reconstruct the phase transitions in the 5'-planes. The 
critical fields Hf''' = 5.7 T, iJf'^''^ = 7.2 T, calculated 
by use of C, and Hgat following the theory y , are shown 
in dashed lines in Fig. [2^. As indicated by the arrow 
diagrams, and in perfect agreement with the diffraction 
experiments, at low fields one expects an incommensu- 
rate spiral structure confined to the (a, b) plane. Be- 
yond the phase transition at i/j^'"^ the spins should form 
a commensurate 3-sublattice configuration with the 2D 
propagation vector (1/3,1/3), as observed experimen- 
tally. At still higher fields, beyond i/g^'"^, the incom- 
mensurate state is expected to be restored. The spins 
will form a modulated fan-type structure, this time os- 
cillating near the field direction in a small angular in- 
terval, still remaining within the (a, 6) plane. Note that 
neutron diffraction failed to detect any incommensurate 
peaks beyond H^^'^'^ . This implies that either the system 
remains disordered, or that the ordering vector is outside 
our search range in reciprocal space and was simply over- 
looked. Besides, the fiop of the spin plane perpendicular 
to the field direction is expected i2i| at iJ = Hgf = 6.5 T. 
Instead of this instability we observe a lost of the mag- 
netic Bragg peaks from S-planes. 



It is remarkable that even as the S-plane go through 
a series of phase transitions, the C-planes remain unaf- 
fected, and vice versa. This implies that magnetic inter- 
actions between each type are direct, rather than medi- 
ated by the intercalated layers of the other type. These 
long-range interactions are likely to be of dipolar origin. 

In summary, we have demonstrated that KFe(Mo04)2 
realizes not just one, but two co-existing and decoupled 
instances of the distorted TLAF model. While one has a 
commensurate ground state and rather simple behavior 
in applied fields, the other features a chiral incommen- 
surate structure and a series of exotic high-field phases. 
Future work should consider the possibility of multifer- 
roic behavior associated with the S-layers. 
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